Supersolid state in fermionic optical lattice systems 
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We study ultracold fermionic atoms trapped in an optical lattice with harmonic confinement by 
combining the real-space dynamical mean-field theory with a two-site impurity solver. By calculating 
the local particle density and the pair potential in the systems with different clusters, we discuss 
the stability of a supersolid state, where an s-wave superfluid coexists with a density-wave state of 
checkerboard pattern. It is clarified that a confining potential plays an essential role in stabilizing 
the supersolid state. The phase diagrams are obtained for several effective particle densities. 



I. INTRODUCTION 



Since the successful realization of Bose-Einstein con- 
densation in a bosonic ^^Rb systenyi, ultracold atomic 
systems have attracted considerable interest One of 
the most active topics in this field is an optical lattice 
system^^i^ii^ which is formed by loading the ultracold 
atoms in a periodic potential. This provides a clean sys- 
tem with quantum parameters which can be tuned in a 
controlled fashion from weak to strong coupling limits. In 
fact, remarkable phenomena have been observed such as 
the phase transition between a Mott insulator and a su- 
perfluid in bosonic systems^. In addition, the superfluid 
stateiS and the Mott insulating stat o^^i^^ have been ob- 
served in the fermionic optical lattices, which stimulates 
theoretical investigations on the quantum states in the 
optical lattice systems. Among them, the possibilty of 
the supersolid state has been discussed as one of the inter- 
esting problems in optical lattice systems. The existence 
of the supersolid state was experimentally suggested in a 
bosonic ^He system^i^ and was theoretically discussed in 
the strongly correlated systems such as bosonic systemsii 
and Bose-Fermi mixtures^^. As for fermionic systems, it 
is known that a density wave (DW) state and an s-wave 
superfluid (SSF) state are degenerate in the half-filled at- 
tractive Hubbard model on the bipartite lattice except for 
one dimension fi^iii which means that the supersolid state 
might be realizable in principle. However, the degener- 
ate ground states are unstable against perturbations. In 
fact, the hole doping immediately drives the system to a 
genuine SSF state. Therefore, it is difficult to realize the 
supersolid state in the homogeneous bulk system. By 
contrast, in the optical lattice, an additional confining 
potential makes the situation different In our previ- 
ous paper,— we studied the attractive Hubbard model 
on square lattice with harmonic potential to clarify that 
the supersolid state is indeed realized at low tempera- 
tures. However, we were not able to systematically deal 
with large clusters to discuss how the supersolid state de- 
pends on the particle density, the system size, etc. This 
might be important for experimental observations of the 
supersolid state in the optical lattice. 

In this paper, we address this problem by combining 
the real-space dynamical mean-field theory (R-DMFT) 



with a two-site impurity solver. We then discuss how sta- 
ble the DW, SSF and supersolid states are in the optical 
lattice system. We also clarify the role of the confining 
potential in stabilizing the supersolid state. 

The paper is organized as follows. In Sec. |TT1 we in- 
troduce the model Hamiltonian and explain the detail of 
R-DMFT and its impurity solver. We demonstrate that 
the supersolid state is indeed realized in a fermionic op- 
tical lattice with attractive interactions in Sec. IIIIl In 
Sec. IIVI we discuss the stability of the supersolid state 
in large clusters. We also examine how the phase dia- 
gram depends on the particle number. A brief summary 
is given in Sec. |Vl 



II. MODEL HAMILTONIAN AND METHOD 

Let us consider ultracold fermionic atoms in the optical 
lattice with confinement, which may be described by the 
following attractive Hubbard model 



H 



it"-*! 



Yv{ri)n^cr, (1) 



where Cia{cl^) annihilates (creates) a fermion at the ith 
site with spin a and = cj^Qcr- t{> 0) is a near- 
est neighbor hopping, U{> 0) an attractive interaction, 
v{r) [= V{r/a)'^] a harmonic potential and the term {ij) 
indicates that the sum is restricted to nearest neighbors. 
ri is a distance measured from the center of the system 
and a is lattice spacing. Here, we define the characteris- 
tic length of the harmonic potential as d = {V/t)~^^^a, 
which satisfies the condition v{d) — t. 

The ground-state properties of the Hubbard model on 
inhomogeneous lattices have theoretically been studied 
by various methods such as the Bogoljubov-de Gcnnes 
equations2^, the Gutzwiller approximation^^, the slave- 
boson mean-field approach^i, variational Monte Carlo 
simulations^^, local density approximation!^ Although 
magnetically ordered and superfluid states are described 
properly in these approaches, it may be difficult to de- 
scribe the coexisting phase like a supersolid state in the 
inhomogeneous system. The density matrix renormal- 
ization group method^S and the quantum Monte Carlo 
mcthod^i are efficient for one-dimensional systems, but 
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it may be difficult to apply them to higher dimensional 
systems with large clusters. We here use R-DMFT— , 
where local particle correlations are taken into account 
precisely. This treatment is formally exact for the ho- 
mogeneous lattice model in infinite dimensions22 and the 
method has successfully been applied to some inhomo- 
geneous correlated systems such as the surface^ or the 
interface of the Mott insulators^!, the repulsive fermionic 
atoma^i^^. Furthermore, it has an advantage in treating 
the SSF state and the DW state on an equal footing in 
the strong coupling regime, which allows us to discuss 



the supersolid state in the optical lattice. 

In R-DMFT. the lattice model is mapped to an effec- 
tive impurity model, where local electron correlations are 
taken into account precisely. The lattice Green function 
is then obtained via self-consistent conditions imposed on 
the impurity problem. When one describes the supcrfluid 
state in the framework of R-DMFT)^ the lattice Green's 
function for the system size L should be represented in 
the Nambu-Gor'kov formalism. It is explicitly given by 
the {2L X 2L) matrix. 



= -tS 



iujncro + {/.J - w {ri)} - % (?w„) 



(2) 



where (Tq, (a = x, y, z) is the ath component of the (2 x 
2) Pauli matrix, ctq the identity matrix, /i the chemical 
potential, u)n = {2n-\-\)'KT the Matsubara frequency, and 
T the temperature. The site-diagonal self-energy at ith 
site is given by the following (2 x 2) matrix. 



(iuJn) 
Si {iUn) 



Si (iw„) 



(3) 



where Si(iix'„) [Si{iujn)] is the normal (anomalous) part 
of the self-energy. In R-DMFT, the self-energy at the 
ith site is obtained by solving the effective impurity 
model, which is explicitly given by the following Ander- 
son Hamiltonian^i^ 

Himp,i = Eikali^^aika + {Dikaik^aikj + h.c.) 

ka k 
ka 

+ ^» X! ^^^<y " 4\ 4i c^i ' (4) 



where o,ikiy{fl\k^) annihilates (creates) a fermion with spin 
a in the effective bath and e,; is the impurity level. We 
have here introduced the effective parameters in the im- 
purity model such as the spectrum of host particles Eik , 
the pair potential Dik and the hybridization Vik- By 
solving the effective impurity model eq. ^ for each 
site, we obtain the site-diagonal self-energy and the lo- 
cal Green's function. The R-DMFT self-consistent loop 
of calculations is iterated under the condition that the 
site-diagonal component of the lattice Green's function 
is equal to the local Green's function obtained from the 



effective impurity model as 



Glat {i^n 



G 



imp, I 



{iuj„) 



When R-DMFT is applied to our inhomogeneous sys- 
tem, it is necessary to solve the effective impurity mod- 
els L times by iteration. Therefore, numerically powerful 
methods such as quantum Monte Carlo simulations, the 



exact diagonalization method, and the numerical rcnor- 
malization group method may not be efficient since they 
require long time to perform R-DMFT calculations. In 
this paper, we use a two-site approximation^Si^, where 
the effective bath is replaced by only one site. In spite 
of this simplicity, it has an advantage in taking into ac- 
count both low- and high-energy properties reasonably 
well within restricted numerical resources^^*^. 

In the two-site approximation, a non-interacting 
Green's function for the impurity model at the ith site is 
simplified as. 



(iuJri) 



ViCF^ 



1 



' iujnCfo — Ei&z — Didx 



where the index k was omitted. The effective parameters 
{Ei, Di,Vi,ei} should be determined self-consistcntly so 
that the obtained results properly reproduce the original 
lattice problem. Here, we use the following equations. 



Re 



V, = 
D, = 



{a~l){n^T^ + Ef + D^) 
b 



1-a' 



(6) 

(7) 
(8) 



where a ~ lm[Q'j' {iluo)]^^ / nT and b = Re[C?"(iLJo)]i2"^. 
Furthermore, the number of particles is fixed in the non- 
interacting Green's function, as 



2T^Re[elp,,(z^„) 



ri=0 



1 

n + 2- 



(9) 



We can determine the effective parameters {Ei,Di, Vi,ei} 
in terms of these equations. 

Here, the effective particle density is defined as p = 
N/ttcP, where N is the total number of particles. This 
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density relates the systems with different sites, num- 
ber of particles, and curvatures of the confining po- 
tentials in the same way as the particle density does 
for periodic systems with different sites. We set t 
as a unit of energy and calculate the density profile 
(riia) ~ '^Tj2n=a^^[^icr («Wn)] + 5 and the distribution 
of the pair potential A.; = 2T^jj^pRe[Fi (iw„)]; where 
Gia{iuJn)[Fi(iu!n)] is thc normal (anomalous) Green's 
function for the ith site. Note that Aj represents the 
order parameter for the SSF state for the ith site. 

In thc following, we consider the attractive Hubbard 
model on square lattice with harmonic confinement as a 
simple model for the supersolid. In this case, it is known 
that the symmetry of the square lattice is not broken in 
the DW, SSF, and supersohd statesi^i^i^. Therefore, 
the point group C^y is useful to deal with the system 
on thc inhomogeneous lattice. For example, when thc 
system with 5513 sites (r < 42.0) is treated, one can 
deal with only 725 inequivalent sites. This allows us to 
discuss the low temperature properties in larger clusters, 
in comparison with those with {d/a = 6.5, N ^ 300) 
treated in our previous paper. 



(a) U/t = 0.0 

(b) U/t = 2.0 



(c) U/t = 3.0 



(d) U/t = 5.0 
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III. LOW TEMPERATURE PROPERTIES 



By means of R-DMFT with the two-site impurity 
solver, we obtain the results for thc system with d/a = 10 
and N ~ 720(/5 2.3). Figures [U and [g] show the profiles 
of thc local density and the pair potential at T/t = 0.05. 
In the non-interacting case {U/t = 0), fermionic atoms 
are smoothly distributed up to r/a ~ 21, as shown in Fig. 
[T](a). Increasing the attractive interaction U, fermions 
tend to gather around the bottom of the harmonic po- 
tential, as seen in Fig. [T] (b). In these cases, the pair po- 
tential is not yet developed, as shown in Figs. [2] (a) and 

(b) , and thereby the normal metallic state with short- 
range pair correlations emerges in the region {U/t < 2). 
Further increase in the interaction U leads to different 
behavior, where the pair potential A^ is induced in the 
region with (riia) ^ 0. Thus, thc SSF state is induced by 
the attractive interaction, which is consistent with the 
results obtained from the Bogoljubov-de Gennes equa- 
tion^. In the case with U/t = 3, another remarkable 
feature is found around the center of the harmonic po- 
tential {r/a < 7), where a checkerboard structure ap- 
pears in the density profile (rii^), as shown in Fig. [1] 

(c) . This implies that the DW state is realized in thc 
region. On thc other hand, the pair potential A^ is not 
suppressed completely even in thc DW region, as shown 
in Fig. [2] (c). This suggests that the DW state coex- 
ists with the SSF state, i.e. a supersolid state appears 
in our optical lattice system. The profile characteris- 
tic of thc supersolid state is clearly seen in the case of 
U/t = 5. Figures md) andll^d) show that the DW state 
of checkerboard structure coexists with the SSF state in 
the doughnut-like region (5 < r/a < 15). By contrast, 
the genuine SSF state appears inside and outside of the 




FIG. 1: (Color online) The density profile (riio-) in the optical 
lattice system with d/a = 10 at T/t = 0.05 when U/t = 
0.0, 2.0, 3.0, 5.0, 7.0 and 10.0 (from the top to the bottom). 



region {r/a < 5, 15 < r/a < 17). Further increase in 
the interaction excludes the DW state out of thc cen- 
ter since fermionic atoms are concentrated around the 
bottom of thc potential for large U. In thc region, two 
particles with opposite spins are strongly coupled by the 
attractive interaction to form a hard-core boson, giving 
rise to thc band insulator with (ui^) ^ 1, instead of the 
SSF state. Therefore, the SSF state survives only in the 
narrow circular region surrounded among the empty and 
fully occupied states. We see such behavior more clearly 
in Figs. [l](f) and[2](f). Note that in the strong couphng 
limit U/t —^ oo, all particles are condensed in the region 
r < rc = d - 1.07d = 10.7a. 

In this section, we have studied the attractive Hubbard 
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(a) U/t = 0.0 
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size, by performing R-DMFT for several clusters with dif- 
ferent d. We here fix J7/t = 5 and ///t ^ —1.58 to obtain 
the profiles of the local particle density and the pair po- 
tential, which are shown in Fig. O We note that the 



(b) U/t = 2.0 
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(c) u/t = 3.0 
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(d) u/t = 5.0 
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FIG. 3: Profiles of particle density (riicr) and pair potential 
Ai as a function of r/d with fixed d — 12.5, 17.5 and 22.5, 
when U/t = 5. 




FIG. 2: (Color online) The pair potential Ai in the optical 
lattice system with d/a = 10 at T/t — 0.05 when U/t = 
0.0, 2.0, 3.0, 5.0, 7.0 and 10.0 (from the top to the bottom). 

model with the harmonic potential to clarify that the 
supersolid state is realized in a certain parameter region. 
However, it is not clear how the supersolid state depends 
on the system size and the number of particles. To make 
this point clear, we deal with large clusters to clarify that 
the supersolid state is indeed realized in the following. 



IV. STABILITY OF THE SUPERSOLID STATE 

In this section, we discuss the stability of the supersolid 
state in fermionic optical lattice systems, which may be 
important for experimental observations. First, we clar- 
ify how low-temperature properties depend on the system 



distance r is normalized by d in the figure. It is found 
that (nicr) and describe smooth curves for r/d ^ 0.6 
and 1.4 ^ r/d ^ 1.7, where the genuine SSF state is re- 
alized. On the other hand, for 0.6 ^ r/d < 1.4, two dis- 
tinct magnitudes appear in {ui^)^ reflecting the fact that 
the DW state with two sublattices is realized. Since the 
pair potential is also finite in the region, the supersolid 
state is realized. In this case, we deal with finite systems, 
and thereby all data are discrete in r. Nevertheless, it is 
found that the obtained results are well scaled by d al- 
though some fluctuations appear due to finite-size effects 
in the small d case. The effective particle density p is 
almost constant in the above cases. Therefore, we con- 
clude that when p ^ 2.3, the supersolid state discussed 
here is stable in the limit with N^d oo. This result 
docs not imply that the supersolid state is realized in 
the homogeneous system with arbitrary fillings. In fact, 
the supersolid state might be realizable only at half fill- 
ing^i^i^i^i^. Therefore, we can say that a confining 
potential is essential to stabilize the supersolid state in 
the optical lattice system. 

Next, we focus on the system with U/t = 5 and 
d/a — 10 to discuss in detail how the supersolid state 
depends on the effective particle density p. The DW 
state is characterized by the checkerboard structure in 
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the density profile {uicr), so that the Fourier transform 
Uq at q = (tt, tt) is appropriate to characterize the exis- 
tence of the DW state. On tlic other hand, the Fourier 
transform Aq at q = (0, 0) may represent the rigidity of 
the SSF state in the system. In Fig. IH we show the 




FIG. 4: (Color online) (n,r7r) / {noo) and Aoo/5(noo) as a func- 
tion of the effective particle density p(= N/nd?) when U — 5t 
and T = 0.05t. A broken line represents the local particle 
density at the center of the lattice in the noninteracting case. 



semilog plots of the parameters normalized by (noo). It 
is found that the normalized parameter Aoo is always fi- 
nite although the increase in the attractive interaction 
monotonically decreases it. This implies that the SSF 
state appears in the system with the arbitrary particle 
number. In contrast to this SSF state, the DW state 
is sensitive to the effective particle density as shown in 
Fig. m These may be explained by the fact that in the 
system without a harmonic confinement (Vq = 0), the 
DW state is realized only at half filling {n = 0.5), while 
the SSF state is always realized. To clarify this, we also 
show the local particle density in the noninteracting case 
at the center of the system as the broken line in Fig. H) 
It is found that when the quantity approaches half-filling 
0.5), {ut^tt) / (noo) takes its maximum value, where the 
SSF state coexists with the DW state. Therefore, we can 
say that the supcrsolid state is stable around this condi- 
tion. Increasing the effective particle density, the band 
insulating states become spread around the center, while 
the DW and SSF states should be realized in a certain 
circular region surrounded among the empty and fully- 
occupied regions. Therefore, the normalized parameters 
(n,r7r) and Aqo are decreased with increase in p. On the 
other hand, in the case with low density p < 0.3, the 
local particle density rii at each site is far from half fill- 
ing even when U/t = 5. Therefore, the DW state does 
not appear in the system, but the genuine SSF state is 
realized. These facts imply that the condition n ~ 0.5 
is still important to stabilize the supersolid state even in 
fermionic systems confined by a harmonic potential.— 

By performing similar calculations for the systems 
with low, intermediate, and high particle densities (p ~ 
0.63, 2.3 and 8.9), we end up with the phase diagrams, as 



shown in Fig. [5] We find that increasing the attractive 
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FIG. 5: (Color online) The phase diagram of the attractive 
Hubbard model on the optical lattice with p ~ 0.63, 2.3 and 
8.9. The density plot represents the profiles of the s-wave pair 
potential as a function of the attractive interaction U/t. The 
DW state is realized in the shaded area. The broken lines 
give a guide to eyes which distinguishes the region with a 
fractional particle density from the empty and fully-occupied 
regions. 

interaction, fermionic particles gradually gather around 
the center of the system, where the empty state is stabi- 
lized away from the center and the band insulating state 
with fully occupied sites is stabilized. It is found that 
the region surrounded among these states strongly de- 
pend on the effective particle density p. The increase in 
the effective particle density shrinks the region, which af- 
fects the stability of the SSF, DW and their coexisting 
states. In particular, the DW region, which is shown as 
the shaded area in Fig. [5l is sensitive to the effective par- 
ticle density, as discussed above. Namely, the local pair 
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potential takes its maximum value around V jt ^ 15, 
which may give a rough guide for the crossover region 
between the BCS-type and the BEC-type states. Wc 
note that the DW state appears only in the ECS region 
iJJ jt ^ 5). This implies that the condition n ^ 0.5 is not 
sufBcient, but necessary to stabilize the supersolid state 
in the attractive Hubbard model with an inhomogcncous 
potential. 

We wish to comment on the conditions to observe the 
supersolid state in the fermionic optical lattice system. 
Needless to say, one of the most important conditions 
is the low temperaturCf^ Second is the tuning of the 
effective particle density p(~ 1), which depends on the 
curvature of the harmonic potential as well as the total 
number of particles. This implies that a confined poten- 
tial play a crucial role in stabilizing the supersolid state 
in fermionic optical lattice systems. In addition to this, 
an appropriate attractive interaction is necessary to sta- 
bihze the DW state in the BCS-type SSF state. When 
these conditions are satisfied, the supersolid state is ex- 
pected to be realized at low temperatures. 

V. SUMMARY 

We have investigated the fermionic attractive Hubbard 
model in the optical lattice with harmonic confinement. 



By combining R-DMFT with a two-site impurity solver, 
we have obtained the rich phase diagram on the square 
lattice, which has a remarkable domain structure includ- 
ing the SSF state in the wide parameter region. By per- 
forming systematic calculations, we have then confirmed 
that the supersolid state, where the SSF state coexists 
with the DW state, is stabilized even in the limit with 
^ cxo, Vb ^ and p ^ const. We have also elucidated 
that a confining potential plays a key role in stabilizing 
the supersolid state. 
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